Valence bond spin liquid state in two-dimensional frustrated spin-1/2 Heisenberg 

antiferromagnets 
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Fermionic valence bond approach in terms of SU(4) representation is proposed to describe the 
Ji — Ji frustrated Heisenberg antiferromagnetic (AF) model on a bipartite square lattice. A uniform 
mean field solution without breaking the translational and rotational symmetries describes a valence 
bond spin liquid state, interpolating the two different AF ordered states in the large Ji and large 
J2 limits, respectively. This novel spin liquid state is gapless with the vanishing density of states 
at the Fermi nodal points. Moreover, a sharp resonance peak in the dynamic structure factor is 
predicted for momenta q = (0, 0) and (jt, 7r) in the strongly frustrated limit J2/J1 ~ 1/2, which can 
be checked by neutron scattering experiment. 
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A spin liquid ground state without any symmetry 
breaking is regarded as one of the most fascinating possi- 
bilities allowed by the physics of two-dimensional quan- 
tum antiferromagnets. Q It is argued that reduced di- 
mensionality, a small spin value, and the presence of com- 
peting interactions may lead to strong enough quantum 
fluctuations to destroy magnetic long-range order (LRO). 
A realistic prototype is the quantum two-dimensional 
spin-1/2 Heisenberg model with the nearest neighbor 
(NN) and next-nearest neighbor (NNN) antiferromag- 
netic (AF) couplings, which has been recently mate- 
rialized experimentally in Li2VOSi04 and Li2VOGe04 
compounds. The model is defined on a square lattice 
by 



H = J % Si • Sj + J2 Si • Sj, 



(1) 



where Ji, J2 > 0, and the AF alignment between spins 
on the NN sites is hindered by the AF coupling of spins 
between NNN sites. 

It is known that, when J2/J1 <C 1/2, the ground 
state has the conventional Neel LRO with magnetic 
wave vector Q = (tt,it), while in the opposite limit 
1 > J2/J1 ^> 1/2, the minimum energy corresponds 
to a collinear state with Q =(±7r, 0) and (0, ±tt), con- 
sisting of two interpenetrating Neel sublattices with in- 
dependent staggered magnetizations. 3 However, when 
J2/J1 ~ 1/2, i.e., in the strongly frustrated limit, the 
degeneracy of the ground state is large, and there is 
a consensus that it corresponds to a spin liquid state 
without LRO. [5j What is the exact nature of this non- 
magnetic ground state turns out to be one of the most 
challenging problems for the physics of frustrated spin 
systems. There have been a number of different pro- 
posals, including the columnal, or spin-Peierls state, 
where the spin rotational symmetry is preserved but the 
translational symmetry is broken, 0, 0,13 the plaque- 
tte state, recovering the x-y symmetry, || a chiral-spin 
state, breaking the PT symmetry. |10| or a truely homo- 



geneous state, not breaking any translational and rota- 
tional symmetries. [Til fl^] For some time the spin-Peierls 
and plaquette states seemed to be favored by numerical 
studi es. |l3tll4lll5| but the most recent numerical simula- 
tions |l6l Il7j on finite size lattices show strong evidence 
against all states breaking translational symmetry, in- 
cluding spin-Peierls and plaquette states, staggered flux 
phase, etc. Moreover, Capriotti et. al, [l|| have ar- 
gued that the corresponding strongly frustrated ground 
state may be characterized by a projective BCS-type 
wave function though the two different AF LRO states 
in large J\ and large J2 limits failed to be reproduced. 
In view of this latest development the earlier proposal of 
homogeneous spin liquid state|lll [12J appears now as a 
promising candidate. How to construct such a state, and 
how can it interpolate the two different AF LRO states 
in two opposite limits is an outstanding issue. 

In this Letter, by using an SU(4) constrained fermion 
representation to describe the spin-1/2 operators on both 
sublattices simultaneously, we propose a fermionic va- 
lence bond (VB) approach to describe the Ji — J2 frus- 
trated Heisenberg AF model on a bipartite square lattice. 
A uniform mean field (MF) solution gives rise to a new 
VB spin liquid state with gapless spin excitations. At 
zero temperature, due to the inter-band nesting between 
two quasiparticle dispersions, the dynamic spin struc- 
ture factor (DSSF) at momenta q =(0,0) and (tt,tt) for 
spins on the same sublattice and on different sublattices 
are found to be equal to each other up to sign, and a 
sharp resonance peak is formed in the strongly frustrated 
limit J%l J\ ~ 0.5. Although the true AF LRO can not 
be reproduced in the weakly frustrated limits, a broad 
peak at (7T,7r) for J%j J\ <C 0.5 is dominant in the static 
spin structure factor (SSSF) of spins on different sublat- 
tices, while for J2/J1 3> 0.5 rather sharp peaks appear at 
(0, ±7r) and (±7r, 0) for spins on the same sublattice. It 
has thus been demonstrated a smooth crossover from the 
Neel to collinear AF (quasi) LRO when changing J2/J1 
from small to large values, and a uniform VB spin liq- 
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uid phase interpolates between the two weakly frustrated 
regimes. 

The Ji — J 2 model on a bipartite square lattice can be 
rewritten as 



P%,1 - + C} A C iA ) + (C/ i2 Ci,2 + Cl 3 Ci^)], 



Bi 



,1 — [\Cl,4pi,l + Ci t xCi t 4,) — (Cl^Cifl + Ci^Ci^)} 



H 



Ji 



+ Jl /] (Sj,A ' S;,A + Sj,B • S;,b) , (2) 
(U) 

where < ij > denotes summation over NN belonging to 
different sublattices, while (i, I) means summation over 
the NN sites of the same sublattice. Each sublattice 
has N sites. As we know, the bipartite lattice struc- 
ture is essential for the spin density wave (SDW) theory 
in describing AF LRO state, where the correspondence 
between the spin operators on the two sublattices is as- 
sumed as S i B — > S i A , S i B — > S ifA , and Sf B — > —S? A , 
reflecting the presence of the Neel LRO. Actually, the 
bipartite lattice structure is also an important setting to 
study various quantum magnetic systems even in the ab- 
sence of AF LRO. However, the correspondence between 
the spin operators on the two sublattices in the VB state 
is no longer the same as that in the SDW theory. 

We generalize the conventional SU(2) constrained 
fermion to an SU(4) representation [18| . The genera- 
tors of the SU(4) fermion representation are given by 
Fg(i) = C} a C hP , satisfying the SU(4) Lie algebra: 
[F$(i),F>?(jj\ = '[8(i,»F?(i) - 8 a , v F$(i)]8 id . The spin 
operators on the sublattice A can be expressed as 



°i,A 



&i,A = ^,2^,1 + ^,4^,3, 



&i,A — 



\cLCi,! - cl 2 a, 2 ) + (c[ 3 c h3 - ct A c iA ) 



while S^ B and S* B are given by interchanging Cj,2 < — > 
Ci,3, reflecting the symmetry of the model. With this 
representation, Sf A and Sf B (a = x,y,z) are proved 
to satisfy their respective commutation relations of the 
SU(2) Lie algebra and [Sf Al Sf B \ = 0. By imposing a 
local constraint ^ p CLC^ = 1 on each lattice site, we 
can further prove that Sf A = Sf B = 3/4. 

In this new representation, the J\ — J 2 model can be 
rewritten as 



Ji \ ■* 



Pi, Pi 



bLb 



i,j i,J 



-qLq 



LI 



(3) 



where the normal ordering has been taken for the first 
and third terms, and four composite VB order operators 
have been introduced 



When uniform VB order parameters are assumed that 
(A id ) = Aic, (B itj ) = -A u , (Pi,i) = A 2c , {Q u ) = 
— A2s, and the local constraint is replaced by a uniform 
Lagrangian multiplier, the MF model Hamiltonian is ob- 
tained as 

H ™f = ^*[{[A-A lc (k)-A 2c (k)]f2 1 

k 

-[A ls (k) + A 2s (k)]ft 2 }* k 

+\ £ $ k {I A - A 2c(k)] «i - A lc (k)n 3 

k 

+A ls (k)n 4 + A 2s (k)f2 2 }$ k 

+2N [ Ji (AL + A 2 ls ) + J 2 (A| c + AL) + A] , 

where the Nambu spinors are introduced 
as * k = (C k!l! C^ 4 ,C_ k> i,C_ M ), 4 = 

(C k 2 , 3 , C_ ki2 , C_k,3), and 4x4 matrices 
arc defined by = a z ® ctq, J~i 2 = a y <g> a x , 



J7. 



CT - 



Moreover, 



A lc (k) = JiAicfcosfc^ + cosfc y ), A 2c (k) = 
4 J 2 A 2c cos fca; cos ky, Ai s (k) = JiAi s (sinfc x + sinfc y ), 
and A 2s (k) = 4J 2 A 2s smk x cosk y . The fermionic 
Matsubara Green function matrices are then derived as 

Gi j4 (k, icj n ) = [iu! n - (A - Ai c (k) - A 2c (k)) Hi 

+ (A ls (k) + A 2s (k))fi 2 ]- 1 , 
G 2 , 3 (k,iw n ) = [iw„-(A-A2c(k))ni + Ai c (k)n 3 

-A ls (k)n4-A2 S (k)n 2 ] _1 , 

whose poles give rise to the quasiparticlc dispersions 



e±(k) = V(A T A lc (k) - A 2c (k)) z + (A ls (k) ± A 2s (k)) 2 

where e+(k) corresponds to the triplet excitations with 
three-fold degeneracy, and e_ (k) corresponds to the sin- 
glet excitations. However, an inter-band nesting prop- 
erty, namely e+(k) = e_(k + Q) with a nesting wave 
vector Q = (jr, tt) is a very important feature for the two 
quasiparticle bands. The usual intra-band nesting for the 
half-filled Hubbard model leads to the AF LRO at zero 
temperature. Similarly, the inter-band nesting will make 
the AF quasi- long range correlations dominant, exclud- 
ing any possible incommensurate density wave states. 

Hence the ground state energy per site is obtained and 
simplified as 

£ = ~jf E MkHJi (AL + Af s )+J 2 (AL + AL)+A. 



The saddle point equations are derived by minimizing the 
ground state energy with respect to Ai c , A 2c , Ai s , A 2s , 
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FIG. 1: The position of the nodal Fermi point (— ko, x , ko tV ) 
in the quasiparticle spectrum e_(k) as a function of J2/J1. 
Due to the inversion symmetry, another nodal point is at 

(&0,X> ky). 



and A. By solving those equations, we can determine the 
saddle point parameters as functions of J2/J1. We notice 
that the two pairing order parameters are competing with 
each other: when J2/J1 < 0.545, Ai s > A 2s , while for 
J2/J1 > 0.545, A ls < A 2s . 

Due to the inter-band nesting property, only the quasi- 
particle spectrum e+(k) is considered for different cou- 
pling parameters J 2 / Ji . Although the uniform VB order 
parameters are assumed above, e+(k) displays two nodal 
Fermi points at (t^o.z, i^o.y), and the quasiparticle den- 
sity of states algebraically vanishes at the Fermi points. 
Moreover, when J2/J1 <C 0.5, the two nodal Fermi points 
are very close to the diagonal line k x + k y = 0, while for 
J2/J1 3> 0.5 the two nodal Fermi points are rotated ver- 
sus the vertical line k x = 0. The position of one nodal 
Fermi point (— ko iX , k , y ) is plotted as a function of J2/J1 
in Fig.l. 

In order to reveal the nature of the new VB spin liq- 
uid state, the DSSFs are calculated. In terms of Nambu 
spinors and <E>j, the spin operators on the two sublat- 
tices are written as Sf A = (^/jfl^^A - < &|r2 5 <I> i )/4 and 

S i,B = (*Jn5*i,A + $ln 5 $i)/4, where fi 5 = a z <g> a z , 
and in the Fourier space the DSSF for spins on the same 
and different sublattices are given by, 



Xls(<3M w m) 



Xi )4 (q,iw m ) ± x 2 , 3 (q^ m ), ( 4 ) 



where 



1 ^ f d*\ 



d 2 k_ 

? 



Tr [n 5 Gi,4(k, iw„)ri 5 Gi,4(q + k, iuj m + iui n )} , 

and X2 zi^i^m) can be obtained by replacing Gi i4 with 
G2,3 in x 14 (q,zw r „). Due to the vanishing density of 
states of quasiparticles around the nodal Fermi points, 
the DSSFs are not divergent at zero temperature, yield- 
ing only algebraically decaying spin-spin correlations 
(quasi-LRO). 

By considering the inter-band nesting, the following 
important relations can be further proved 

XA,A(Q, iuJ m) = XlsCQ^m) 



FIG. 2: The imaginary part of the dynamic spin structure 
factor 7mx(Q,w > 0) for different values of J2/J1. 



= XA,A(0,«W m ) = -XA,B(°, iuJ m) = X{Q,iUm), (5) 

giving rise to the characteristic features of frustrated AF 
correlations of spins on the bipartite square lattice! After 
analytic continuation, its imaginary part of the correla- 
tion spectrum at T — is approximated to be 

Im X (Q,c > 0) « ^ J 0^5 [u C+ (k) - C _(k)] 

f [A - A 2e (k)] 2 - Af c (k) Af s (k) + AL(k) 1 
1 e + (k)e_(k) /' 

which is delineated in Fig. 2. Increasing the coupling pa- 
rameter J 2 /Ji, a sharp resonance is gradually formed in 
the strongly frustrated limit J 2 /Ji ~ 0.5 and disappears 
quickly away from it. The resonance peak simultane- 
ously appears at momenta (7r,7r) and (0,0), and may be 
regarded as a signature of the present fermionic VB spin 
liquid state in the strongly frustrated limit of the model, 
which is ready to be verified by future experiments. 

Moreover, to compare with the AF LRO states in the 
large Ji and large J 2 regimes, respectively, the SSSFs at 
T = are evaluated. A smooth crossover is obtained 
for the SSSFs XaaM an d Xab^) when increasing the 
coupling parameter through the strongly frustrated limit 
J2/J1 = 0.5. In Fig. 3, for J 2 /Ji = 0.3 a rather broad 
peak is clearly identified at q = (it, n) in \\ s( t l)' while 
Xa a (l) has no distinctive features. This peak manifests 
the dominant quasi-long range AF correlations towards 
the Neel LRO in the limit of small J2 / J\ ■ In contrary, in 
Fig.4 for J2/J1 = 0.7, the SSSF Xa.aM displays rather 
sharp peaks at momenta (±7r, 0) and (0, ±7r), while the 
SSSF Xa b (l) I s very small and spreads over a certain 
range. These peculiar features indicate the quasi-long 
range AF correlations towards the collinear LRO in the 
limit of large J2I J\- Therefore, it has been shown that 
a fermionic VB spin liquid state in terms of the SU(4) 
representation exists and approximately interpolates be- 
tween the two AF LRO states in the large J\ and large 
J 2 limits, respectively. 

In conclusion, an SU(4) constrained fermion represen- 
tation has been used to describe the J\ — J 2 frustrated 
Heiscnberg AF model on a bipartite square lattice, and 
a fermionic uniform MF solution gives rise to a novel VB 



4 




FIG. 3: The static spin structure factor at zero temperature 
Xa,b(cl) f° r J2/J1 — 0.3. A broad peak at (tv,tt) displays 
quasi-long range AF correlations in the large J\ limit. 



spin liquid state with gapless spin excitations, not break- 
ing any translational and rotational symmetries. This 
MF solution has reproduced correct asymptotic behavior 
in the weakly frustrated limits with the provision that the 
true LRO is approximated by quasi-LRO in our scheme. 
To the best of our knowledge, no other theoretical ap- 
proach has succeeded in doing so in this complicated 
problem. The success of our MF theory is mainly due 
to the appropriate choice of representation which grasps 
the essential features of the strongly frustrated state. Of 
course, the eventual "proof" of our theory should come 
from the experiment. In particular, our explicit predic- 
tion of the sharp resonance peak at momenta (tt, tt) and 
(0, 0) in the imaginary part of the DSSF in the strongly 
frustrated limit J2/J1 ~ 0.5 can serve as a crucial ex- 
perimental test. There are still many open problems, for 
example, what is the nature of transitions between this 
VB state and AF LRO states, their possible coexistence, 
etc. These issues require further studies. 




FIG. 4: The static spin structure factor at zero temperature 
Xa a(°l) f° r J2/J1 = 0.7. Sharp peaks appearing at (±7r, 0) 
and (0, ±7r) indicate the strong quasi-long range correlations 
of the collinear AF state. 
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